function [ y ] = bvalue( f_D_pD_est, f_D_pD_tr, b_as)
%BVALUE Summary of this function goes here
%   Detailed explanation goes here
y=zeros(1,length(b_as));

for i=1:length(b_as)
    b=b_as(i);
    tmp_f=@(f_D_pD_tr) (1-f_D_pD_tr(1))*exp((-b)*f_D_pD_tr(2)) + f_D_pD_tr(1)*exp((-b)*f_D_pD_tr(3));
    if floor(max(b_as(i))/100)-2>0
        f_D_pD_tr(1)=f_D_pD_tr(1)+0.005;
    end
    y(i)=tmp_f(f_D_pD_tr);
end
 plot(b_as,log(y),'r');
% y=tmp_f(f_D_pD_tr)+(rand()*(c+c)-c)*tmp_f(f_D_pD_tr);
% f=@(f_D_pD_est) ((1-f_D_pD_est(1))*exp((-b)*f_D_pD_est(2)) + f_D_pD_est(1)*exp((-b)*f_D_pD_est(3))-y)';
% [f_D_pD]= LMFnlsq2(f,f_D_pD_est,'Display',[0,0]);

% plot(b,log(y),'o',b,log(tmp_f(f_D_pD)),'r');
% plot(b,log(tmp_f(f_D_pD_est)),'g',b,log(tmp_f(f_D_pD_tr)),'r');
end

